Are back-to-back particle— antiparticle correlations observable 
in high energy nuclear collisions? 
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Analytical formulae are presented which provide quantitative estimates for the suppression of the 
anticipated back-to-back particle-antiparticle correlations in high energy nuclear collisions due to 
the finite duration of the transition dynamics. They show that it is unlikely to observ the effect. 
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I. INTRODUCTION 

In 1996 Asakawa and Csorgo [JJ suggested that back- 
to-back particle-antiparticle correlations should be ob- 
servable in high-energy nuclear collisions, see also earlier 
considerations in Refs. [2,3. Based on a sudden transi- 
tion assumption huge effects were predicted for this phe- 
nomenon. About a further dozen applications followed. 
In a subsequent paper together with Gyulassy [4] it was 
shown that the finite duration r of the transition re- 
duces the effect. In that paper, however, the authors 
used a discontinuous and therefore unrealistic transition 
profile that led to large ultraviolet Fourier components 
and thus to a very moderate reduction of the effect of 
order (2u>t)~ 2 . Here 2uj is the energy required to pro- 
duce the pair. In this short note it will be shown that 
with a smooth transition profile and realistic transition 
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FIG. 1: Optical ray pictures in space-time. Left: the standard 
reflection-transmission situation, where a wave traverses from 
one medium (yellow) to another one (e.g. vacuum) at a spa- 
tial interface. Right: the sudden transition in time from the 
medium to the vacuum; here the frequency ui is discontinuous, 
creating a second component with negative oj which has the 
interpretation of an antiparticle with opposite momentum. 
The (blue) dashed line separates the media, the arrow sense 
distinguishes between particle and antiparticle. By means of 
Lorentz transformations the concept can directly be general- 
ized to hyper-planes and even to curved hyper-surfaces. 
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durations the suppression is rather exponential, oc e 4wT , 
deferring to observe the effect. 



II. THE SUDDEN PICTURE 

The main assumption of the original approach [JJ is 
the sudden change from the in-medium situation to that 
in vacuum. The authors used the Bogolioubov-Valatin 
(BV) transformation to describe the effect, a picture 
that may not be so intuitive for many of us. Indeed 
the effect is nothing else than a standard "reflection- 
transmission" problem at the interface of two media, a 
well known problem in physics, in particular in optics. 
Thus, transcribing our wisdom from the spatial situa- 
tion, like a wave traversing from one medium to another 
one at a sharp interface (see Fig. IT] left) to the here con- 
sidered sudden transition in time (Fig. [lj right), precisely 
recovers the results presented in [TJ The common 
boundary condition is that one has one incoming wave 
and two time-forward propagating outgoing wave com- 
ponents. As in Refs. [TJ 0] we restrict the discussion 
to (charged) rclativistic bosons described by the Klein- 
Gordon (KG) equation. For the sudden case the wave 
functions with positive and negative energies prior and 
post to the sudden transition can then be written as (us- 
ing units with h=c=l) 
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Here Jlk and Wk are the single-particle energies at mo- 
mentum k in the two media. Continuity of vf'(x) and of 
d^(i)/dt at transition time t = determines the coeffi- 
cients on the vacuum side to 
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while due to the spatial homogeneity the spatial momen- 
tum k remains unchanged. Within the standard rela- 
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tivistic quantum field theory (RQFT) interpretation of 
the negative energy components as antiparticles with op- 
posite spatial momenta, cf. Fig. [T] right frame, this result 
identically reproduces that given in Refs. [TJ H] within 
the BV picture. In particular the ratio of the antiparti- 
cle over particle component on the vacuum side 
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is the analog of the well known reflection coefficient for 
a wave traversing from one medium to another (Fig. [TJ 
left) 
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where K and k denote the moduli of the spatial mo- 
menta in the two media. In the spatial case solely kj_ is 
discontinuous (w and k|| are continuous). In the sudden 
case, however, k is continuous, while the discontinuity in 
time causes a discontinuity in ui thereby creating a sec- 
ond component with negative ui, i.e. an antiparticle with 
opposite spatial momentum k (Fig. [TJ right frame). This 
mechanism led to the back-to-back particle-antiparticle 
correlation picture advocated in Ref. [TJ. Besides small 
antiparticle components arising form the existing parti- 
cles in the medium, the time dependent interaction can 
also create genuine particle-antiparticle pairs out of the 
vacuum. The latter, which is the bosonic analog to the 
Dirac case, where a time-dependent interaction pulse can 
lift a particle from the filled Dirac sea to the particle 
space 1 , can formally be included by adding a to 
the boson occupations n k of the antiparticles. After the 
transition the one-body density becomes 

m(k) = |c k | 2 n k + |s k | 2 (7i k + l) (8) 

" k = exp(|fi k |/T) + l' (9) 

while the back-to-back correlation function of particle- 
antiparticle pairs over the product of single yields be- 
comes [JJ 
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Due to the sudden pair creation processes this back-to- 
back correlation ratio can attain huge values, once the 



1 In the Dirac case also spatial transitions (Fig. [TJ left) lead to 
spontaneous pair production, known as Klein paradox [5], once 
the vector potential changes by more than twice the rest mass. 



statistical occupations ro k fall below |s k | 2 . The subse- 
quent considerations will show that this effect is an arti- 
fact resulting from the sudden limit. 



III. THE CONTINUOUS TRANSITION CASE 

Compared to the BV transformation considered in [TJ 
E] the wave dynamical picture used here can easily be 
generalized to the continuous transition case by solving 
the time dependent Klein-Gordon equation (suppressing 
in the following the dependence on spatial momentum k) 

d 2 V(t) +n 2 (t)V(t) = (13) 

for a smoothly time-dependent dispersion relation 

n 2 (t) = lo 2 + n R (i) =uj 2 + n R (-oo)F(t). (14) 

Here negative times refer to the situation in the medium 
where II R denotes the time-dependent retarded polariza- 
tion function 2 , while the vacuum case is attained for large 
positive times. Thereby F(t), cf. Fig. [2] below, will later 
be used to parameterize the transition profile. 

Besides direct numerical solutions of ( fl3| , there are 
approximate analytical methods to obtain the reflection 
coefficient. These concern the single reflection limit as 
well as refined semi-classical methods. 

A. Single reflection approximation 

A simple generalisation of the sudden limit towards a 
contineous treatment is provided by the single reflection 
approximation (SR). In this scheme the continuous func- 
tion F(t) is approximated by a sequence of steps each 
providing a reflected wave component according to (JsJ) . 
In the continuum limit the coherent sum of these partial 
reflections leads to 
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dt m) exp (r 2i l! dt ' m ) (15) 

dt mt)Q(t)^i(t). (16) 



Thereby the appropriate phase coherence is approxi- 
mately provided by describing the unreflected wave com- 
ponents of particles (i) and their antiparticle partners (f ) 
in the semi-classical (WKB) limit as 



*i(t) 



exp(-ij^dt'n(t'))/y/2n{t). (17) 



2 In general the polarization function II can also be non-local in 
time and thus be energy dependent in the semi-classical interpre- 
tation. Then the mapping from the medium to the vacuum case 
will involve corresponding "z"-factors, z = 1/(1 — c9Re U 11 /du) 2 ), 
as e.g. shown in Kadanoff-Baym final-state interaction picture 
[6]. For the example cases discussed here we will discard such 
generalization and stick to time local cases. 
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Since f2 is peaked and limited to a narrow range in t, the 
antiparticle amplitude A is essentially given by the time 
Fourier transformed of Cl(t), i.e. 
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with abbreviations f2 = f2(— oo) and a; = f2(+oo). There- 
by w = f2(t) is chosen around the maximum of f2. 
As multiple reflections at the different steps are sup- 
pressed, the approximation solely recovers terms linear 
in Afl = Q — ui. Compared to the sudden result (6j) the 
suppression of the antiparticle amplitude caused by the 
finite duration of the transition is then given by g(2Q). 
This was also the approximation scheme used in Ref. [3] 
where, however, a completely unrealistic form for f(t) 
was considered, namely f(t) = 8(t) exp[— </t]/t. This 
form contains a sharp jump which leads to unrealisti- 
cally high Fourier components for large w, and thus only 
to a very moderate suppression of the effect! More re- 
alistic forms for Q(t) generally will lead to exponential 
suppression factors, as e.g. for the following analytically 
solvable cases 
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(20) 

In order to assure comparable results all functions are 
chosen such that F(t) monotonically falls from 1 to 
with a maximum time derivative at t = of /(0) = 
-F(0) = l/(2r), cf. Fig. ^ Then r is approximately the 
"half time" for the change in f2 2 (t). For short transitions 
times r — > 0, one verifies the sudden result ([6| since then 
g(2uj) — > 1. For large wt all three cases indeed lead to 
exponential suppression factors 
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FIG. 2: (Color online) -F(t) for the three example cases of 
table <\20\, (a) red dashed, (b) blue dotted, (c) black full line. 



with C = (1,4, (47twt) 2 ) and a = (2/tt, 1, tt/2) for the 
cases (a) to (c) of (20|. Thereby the asymptotic expo- 
nential behavior is directly determined by the imaginary 
part of the nearest complex pole position of F(t). 

The physical case where the polarization function is 
switched on and off, i.e where F(t) takes the form of 
a bell shape can simply be retrieved by replacing F(t) 
by its time derivative F^ e \\(t) = 2rf(t) with f(t) from 
table ( 20 ) . The corresponding "pulse" has a maximum 



of 1 and a time integral of 2r. The suppression factor 
( |21[ ) for the pair production then will simply attain an 
additional pre-exponential factor of (4cDt) 2 . 



B. Complex semi-classical method (CWKB) 

An alternative method to access the reflected wave 
component is provided by the complex WKB method [7] . 
In that case the "reflected" wave component results from 
the analytic continuation of the "unreflected" WKB wave 
^i, cf. Eq. (17), around the nearest complex valued turn- 

0. 



ing point to where f2 2 (£o) = 0. Then the reflection coef- 
ficients 1Z or A are given by [7] 

A CWKB = exp[-2|Im S(t ,t)\] (for real t), (22) 

where 



S(t ,t) = IdHy/mf), 



(23) 



is the action integral from the nearest complex turning 
point to to some time t on the real axis. Since the action 
integral along the real axis does not contribute to Im S, 
a rough estimate can simply be obtained by choosing the 
integration contour starting from to just parallel to the 
imaginary axis till the real axis. Assuming that along this 
contour the integrand can essentially be approximated by 
its value on the real axis one obtains the estimate 



ImS(to) « ^(Reto) Im t 



(24) 



For the three analytic cases of table ( 20 ) one verifies that 
the complex turning points are located close to the pole 
positions of F(t) with values for Im to = ar with a — 
2/tt , 1 and 7r/2, respectively for the cases (a) to (c) in 
( 20 1 . Thus the pair creation rate then becomes 



(25) 



in agreement with the leading exponential terms of the 
single reflection approximation given in (21). 



The exact Fermi function case 



The Fermi- function case ( 20 ;) is particularly interest- 
ing, since the corresponding KG equation ( 13 ) can be 



solved in closed form, providing the exact result [8 
sinh(7r(fi - w)r/2) ° 
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(26) 
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It generalizes the sudden result (|6| to the continuous 
Fermi-function case ( 20 ;) and further confirms the limit- 



ing cases of the single-reflection approximation ( 18 ) and 



the CWKB result ( 25 ) within their validity ranges. Also 



the corresponding action integral ( 23 1 and thus its imag- 



inary part can be obtained in closed form 

■SVcrmi^Oj t) = Qt art anh - — lot artanh — — (27) 

i I LU 

|Im S*Fcrmi(to,t)| = \ Min(n, w) t (for real t), (28) 

the imaginary part resulting from that artanh function 
with argument larger than 1. The corresponding CWKB 
amplitude ( 22 ) then becomes 



CWKB I 2 

Fermi I 



= exp [-2-7T Min(0, w) t] 



(29) 



in full compliance with the semi-classical limit {lot 3> 1, 
ASlr 3> I) of the exact Fermi- function result (26). 



IV. 



SUMMARY AND CONCLUDING 
REMARKS 



For nuclear collisions with a freeze-out temperature T 
and a correspondingly smooth time dependence of the 
polarization function n, the back-to-back correlation ( 12 ) 
can then be estimated to 

|c 2 (k, -k) - 1| « 0(|s k | 2 A*k) « e - n ^ 4aT - 2 / T \ (30) 

valid for |sk| 2 <§C «k <C 1. Thus, the original effect 
advocated in Refs. [U 0] of unlimitedly large correla- 
tions with increasing resulting in the sudden case 
(t = 0) turns into the opposite behavior, once the du- 
ration time t exceeds the inverse temperature 1/T. Al- 
ready for charged pion pairs (ir + ir~) with fik ~ 300 MeV 
the correlation (30) falls below the 10 -5 level for tran- 



sition times t beyond 3 fm/c, at typical freeze-out tem- 
peratures of T ~ 140 MeV. Since optical potentials are 
generally proportional to the density, the characteristic 
transition-time scale is rather given by the expansion- 
time scale [SHU] in nuclear collisions. Its full-width-half- 
maximum values generally lie beyond 4 fm/c, leading to 



even larger suppression factors. Thus, the pair creation 
process due to the time variation of the mean field be- 
comes unmeasurably small in nuclear collisions. 

The discussed particle-antiparticle correlation effect 
rests on a coherent single-particle picture. Among oth- 
ers it largely ignores collisional effects that are known to 
dominate the nuclear collision dynamics. Thus, in in- 
dividual events the one-body field will depart from the 
ensemble mean due to fluctuations caused by stochastic 
processes. Such microscopic processes can involve much 
shorter time scales which, however, are also ultra vio- 
let restricted by the frequencies accessible in the system. 
For thermal systems this limit is given by the inverse tem- 
perature scale, i.e. r > 1/T, limiting the "thermal" pro- 
duction of hard probes, such as the particle-antiparticle 
pairs. Such stochastic processes generally add incoher- 
ently to the two-particle yield. Therefore they can ap- 
propriately be included in transport codes by simulating 
the corresponding microscopic collision processes. The 
resulting pair correlations, however, are no longer strictly 
back-to-back, as they arc influenced by thermal motion. 

The question, whether a process can be treated in the 
sudden limit or not, depends on its intrinsic quantum 
time scale resulting from the uncertainty principle to 
tq s» 1/AE, where AE is the energy transfer. Driven 
by a certain field, the dynamics is "sudden" , if the 
typical time duration r, during which the field changes, 
is short compared to tq: then the wave functions stay 
"inert" across the transition thereby creating several 
components in the eigenstates of the "new" Hamiltonian. 
In the opposite limit t > tq the quantum states change 
and adjust adiabatically. As shown here for the pair 
creation case, where AE = 251k, the creation of states 
involving an energy transfer AE 3> 1/r then becomes 
exponentially suppressed. 
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